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Abstract 

The elastic interaction, induced by the lattice distortion due to the difference of the molecular 
size, causes an effective long-range interaction. In spin-crossover (SC) compounds, local bistable 
states, i.e., high-spin and low-spin states have different molecular sizes, and the elastic interaction 
is important. In bipartite lattices, e.g., the square lattice, the ground state can be two types 
of phases: ferromagnetic-like and antiferromagnetic-like phases. In systems like SC compounds, 
the former phase consists of all small or large molecules, and the latter phase has the configu¬ 
ration of alternating small and large molecules. In fact, both cases are observed in SC systems. 
In this paper, we have studied the effect of the short-range interaction in the elastic system on 
the properties of those order-disorder phase transitions. We have obtained a phase diagram in 
the coordinates of the temperature and the strength of the short-range interaction, including the 
metastable structures. We show that effects of the short-range interaction are essentially differ¬ 
ent for ferromagnetic-like and antiferromagnetic-like phase transitions. In the ferromagnetic-like 
transition, the long-range interaction of elasticity is relevant, and the system exhibits a phase 
transition in the mean-filed universality class. In this case, the long-range interaction strongly 
enhances the ferromagnetic-like order, and it works cooperatively with the short-range interaction. 
In contrast, in the antiferromagnetic-like transition, the elastic interaction slightly enhances the 
antiferromagnetic-like order, but essentially it does not contribute to the ordering, and the system 
shows a transition in the Ising universality class. We have found that in the border region between 
ferromagnetic-like and antiferromagnetic-like phases, the antiferromagnetic-like phase has an ad¬ 
vantage at finite temperatures. We discuss the critical properties of two-step SC transitions with 
comparison between the elastic interaction model and conventional SC models (Ising-like models). 

PACS numbers: 75.30.Wx 68.35.Rh 64.60.De 75.60.-d 
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I. INTRODUCTION 


Spin crossover (SC) compounds have attracted much attention in a wide variety of phase 
transitions as well as their potential applications.^^^The spin-crossover system has bistable 
states, i.e., a low spin (LS) state and a high (HS) state with an entropy difference between 
the two states. Entropy-induced phase transition and other cooperative phenomena have 
frequently been studied by the Ising model with the different degeneracies between the LS 
and HS states (called Ising-like model) although the Ising interaction was only introduced 
as the simplest description of the cooperative interaction.!!^ From the microscopic viewpoint, 
the size of each molecule in SC solids changes depending on the spin state, LS or HS, and the 
importance of the elastic interaction has been suggested for the cooperative interaction.!!^® 
By treating the change of the molecular size explicitly, it was shown that the elastic inter¬ 
action, caused by the lattice distortion due to the difference of the molecular sizes between 
the LS and HS states, leads to the cooperative interaction for the SC phase transition,® 
where the effective long-range interaction of elastic origin induces a phase transition. In 
this limiting case without any short-range interaction, the nature of the phase transition has 
been studied from the view point of the pure long-range interacting model, and it was clari¬ 
fied that the phase transition belongs to the mean-field universality class.® Other important 
features of the transition have also been clarified with the use of this kind of modeling®®. 

In realistic compounds, the short-range interaction also plays a role in phase transitions. 
The potential energy function between molecules is considered to depend on the spin states 
of the molecules.®® Then, this dependence is expressed by using a kind of short-range 
interaction of the spin states. In this case, the interaction between molecules consists of 
both short-range and long-range components. 

The importance of the short-range antiferromagnetic-like interaction has also been 
suggested in SC phenomena in the context of explaining two-step transitions.®® An 
antiferromagnetic-like phase, in which the LS and HS molecules align alternately, is re¬ 
alized as a thermodynamic phase in the middle temperature regions between the LS phase 
at lower temperatures and the HS phase at higher temperatures. 

In bi-nuclear systems, i.e., two sites in a unit molecule, the alternate structure may be eas- 
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ily understood. However, the alternate structure appears also in mono-nuclear compounds, 
although the crystal has no sublattice symmetry. In theoretical studies of two-step SC 
transitions, a short-range antiferromagnetic-like interaction was adopted into conventional 
models such as Ising-like models.®®^. (In Sec. an outline of two-step SC transitions and 
a related discussion are given.) 

Competition and interplay between the short-range and long-range interactions are in¬ 
teresting topics in phase transitions.^^® In the pure short-range model, clustering of the 
ordered phase takes place near the critical temperature, leading to the divergence of the cor¬ 
relation length of the order parameter. However, the long-range interaction suppresses the 
generation of domains, and the conhguration is uniform even at the critical temperature.^ 
The crossover from the short-range Ising universality class to the mean-field universality 
class was studied for the case of a ferromagnetic short-range interaction.^^SMl those works, 
it was found that the long-range interaction is always relevant, and the system exhibits a 
phase transition in the mean-held universality class. 

Elastic interactions with a size difference between atoms have been also studied for al¬ 
loy system^Si®^, where elastic potential energies with empirical parameters were adopted for 
neighboring atomic species (Si-Si, Si-Ge, Ge-Ge), and a bond angle potential was introduced 
to maintain the diamond structure. In those models different elastic constants for different 
atomic species lead to a kind of short-range interactions. The critical properties were found 
to be mean-held like for ferromagnetic-like order, while they were suggested to be of rigid 
Ising type for antiferromagnetic-like order, although the transition temperature is signih- 
cantly diherent from that of the Ising model. However, for the ferromagnetic-like order, the 
elastic interaction of these models causes asymmetry of the entropy ehect between Si-Si and 
Ge-Ge phases due to diherent elastic constants for diherent atomic species. Asymmetry of 
the entropy ehect between two phases is also induced by the bond angle potential as we 
discuss in Appendix Indeed, in those studies an artihcial held was applied to avoid the 
asymmetry and to constrain the systems to the coexistence line. 

In the present work, we use the same elastic constant for diherent species. To maintain 
the lattice structure (square lattice), we adopt next nearest neighbor interactions, in which 
the asymmetry of the entropy ehect is negligible as we see in Appendix We focus on 
the dependence of the critical phenomena on the parameters of the short-range interaction. 
Our modeling enables us to study systematically the ehect of the short-range interaction 
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on the critical phenomena without applying an additional held to constrain the systems 
to the coexistence line. Here we can avoid not only the ambiguity and complexity of the 
asymmetry of the entropy effect but also contamination of the short-range interactions due 
to different elastic constants for different molecular species. 

In this paper, we systematically investigate the effects of both ferromagnetic-like and 
antiferromagnetic-like short-range interactions on the ordering process in a unihed model 
and present a phase diagram in the coordinates of the temperature and the strength of the 
short-range interaction, including the metastable structures. We hnd that the nature of the 
critical phenomena is different in the cases of ferromagnetic-like and antiferromagnetic-like 
transitions. We clarify that the contributions of the short-range interaction and the elastic 
interaction to enhancement of the ordering are essentially different in the ferromagnetic-like 
and antiferromagnetic-like transitions. We investigate in detail the border region between 
ferromagnetic-like and antiferromagnetic-like phases and hnd that the antiferromagnetic-like 
phase has an advantage at nonzero temperatures. To expand also discussion of the critical 
properties of two-step SC transitions, we briehy summarize the characteristic features of the 
types of SC transitions from the view point of the phase diagrams of the Ising-like models 
and examine the difference of the properties between the Ising-like models and the elastic 
interaction model. We also analyze several kinds of interactions to maintain the lattice 
structure (symmetry), focusing on asymmetric entropy effects for broken symmetry phases 
due to the difference of molecular sizes. 

The rest of the paper is organized as follows. In Sec. the model and method 
are presented. In Sec. |III[ we discuss the critical properties of ferromagnetic-like and 


antiferromagnetic-like order parameters. In Sec. IV, we show the phase diagram of the 
present model in the coordinates of temperature T vs. the short-range interaction param¬ 
eter J. In Sec. |V| we show an outline of two-step SC transitions. In Sec. VI, we discuss 
the critical properties of two-step SC transitions and give summary and discussion. In Ap¬ 
pendix 1^ we give a discussion about types of interactions to maintain the square lattice 
and about the symmetry between HS and LS. 
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FIG. 1: (color online). ( a) Configuration of the ferromagnetic-like phase (low-spin phase), 

(b) Conhguration of the ferromagnetic-like phase (high-spin phase), (c) Configuration of the 
antiferromagnetic-like phase. Blue and red denote low-spin and high-spin states, respectively. 
High-spin molecules are larger than low-spin molecules. 

II. MODEL AND METHOD 

We consider a system which consists of molecules on a square lattice. The positions of 
the molecules {vt} and also their radii can be changed, and then the lattice can be distorted. 
We assume that the square-lattice topology is not broken. We adopt the following elastic 
interactions between the nearest neighbor molecules, 


^nn = y Ki “ 


( 1 ) 


(»d) 


where r*j is the distance between the centers of the Hh and jth molecules. Each molecule 
takes the low spin (LS) state (cTj = —1) or the high spin (HS) state (cTj = 1) and its radius R 


depends on the state. The LS molecule has a smaller radius than the HS molecule; < Rh 


where i?L (-Rh) is the radius of the LS (HS) molecule. When the molecules contact each 
other, the energy of R^n is minimum. In Fig. [T| the configurations of the minimum energy 
are depicted. 

To avoid a global deformation to rhombic shape and also to maintain the square lattice, 
a small perturbation, such as the following next nearest neighbor interaction, is necessary. 



( 2 ) 


((i,k)) 

Then, the system of the elastic interaction has the form 



( 3 ) 
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We set k 2 = fci/lO in this study. 

Now, as the short-range interaction, we introduce the following nearest-neighbor Ising 
interaction 

Uis = (4) 

(*j> 

Taking into account the energy difference D and the ratio g of the degeneracies between 
the HS and LS states,®Eil we also add the following term (see also Sec. 0 . 

"HtS = i (£>-Coring) (5) 

i 

Thus, the Hamiltonian of the system is given by 

"^tot = "^Elastic + HlS + Heff- (6) 

In the present study, we focus on the critical phenomena. Thus, we need to investigate 
the system along the coexistence line of the model. It is approximately realized if we set 
D — k^T\i\g = 0, i.e., Hefr = 0 (equivalent to D = 0 and g = 1). We give an analysis about 
the coexistence in Appendix]^ We set the parameter values as = IT, = 1, =40 

and /c 2 = 4. 

Here we apply a Monte Carlo (MC) method with NPT ensembl^^, where the pressure 
is set P = 0, for the square lattice (2D) with periodic boundary conditions. In the MC 
method, we choose a molecule at site i, and update the spin state ai and the position of the 
molecule {xi,yi). Then we update the volume of the total system under P = 0. One Monte 
Carlo step (MCS) is dehned as L x L times of these procedures, where L denotes the linear 
dimension (number of sites) of the system. 

A. Ground state configuration 

As mentioned above, some additional interaction is necessary to maintain the square 
lattice. We here adopted the form (Eq. ([^ ). Other types of interactions can also be used 
to maintain the shape (see Appendix [A]). 

Within Pnn, the ferromagnetic-like conhguration (Fig. [^a) and (b)) and the 
antiferromagnetic-like conhguration (Fig. [^c)) are degenerate. However, "Hnnii resolves this 
degeneracy. The total energy of the ferromagnetic-like state per molecule is = 
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and that of the antiferromagnetic-like state is 

^ = Y + (7) 

where is the coordination number (z = 4). The energy difference between the two states 

is 

AE ,_k^(R„-RLfz 

L2 - 12 -■J^+ 2 ' ' ' 

Therefore, for J = 0, AE > 0 and the ferromagnetic-like state is favored. Substituting 
Rh = 1.1, Rl = 1-0, z = 4, and /c 2 = 4, AE/ = 4{J + 0.02). Thus, at J = Jo = —0.02, 
the ground state changes between the ferromagnetic-like (J > Jo) and antiferromagnetic-like 
(J < Jo) states. Hereafter, we define Jo as the origin of J. 


III. CRITICAL PROPERTIES OF TWO ORDER PARAMETERS 


We study the dependence of the critical properties of the model on the short-range inter¬ 
action J. In the present model, the magnetization (m) and staggered magnetization (mgg) 
are the essential order parameters. The dehnitions of m and rrigg are given by 




L2 


(9) 


and 


rrisg = 


E,(-i) 


ix-\-ix 


(Ti 


L2 


( 10 ) 


where {ix, iy) is the integer coordinate of the ith molecule which numbers the 2D lattice. It 
should be noted that m {msg) is not real (staggered) magnetization but a kind of pseudo 
(staggered) magnetization to show ferromagnetic-like (antiferromagnetic-like) order. 

In order to study the critical phenomena, we analyze Binder cumulants for both order 
parameters. Binder cumulants for ferromagnetic-like and antiferromagnetic-like orders are, 
respectively, dehned as 




{m^)L 

3(m2)2 


and U^^{L) = 1 — 


{mtg)L 

3(^sg)i‘ 


( 11 ) 


At the critical temperature, the Binder cumulants for different values of L cross, and the 
value at the point depends on the type of phase transition. 



We also investigate the correlation function of the spin state in the vicinity of the 
critical point to catch the difference of the ordering between ferromagnetic-like and 
antiferromagnetic-like phases. The dehnition of the correlation function is given by 




( 12 ) 


Here ly) denote the position of the /-th spin, and i and j are taken in the range 0 < 



FIG. 2: (color online) (a) A snapshot of the configuration at T = 0.24 (near = 0.238) for 
J = 0.01 and L = 50. A uniform configuration is seen, (b) Temperature dependence of Uj 
for various system sizes L when J = 0.01. The horizontal dotted line {Uj = 0.27) denotes the 
fixed-point value of the Binder cumulant for the mean-field universality class. 


A. Ferromagnetic-like phase transition 

First we study the case of the ferromagnetic-like phase transition. Here we consumed 
1000,000-2000,000 MCS for the initial equilibration and the following 1000,000-6000,000 
MCS at each temperature to obtain physical quantities. 

Figure (a) shows a typical conhguration near the critical temperature T^. when J = 
0.01. The red (blue) solid circles denote HS (LS) molecules. The characteristic features 
are the same as we studied for J = 0,^^ and no clustering occurs. In Fig. (b), the 
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Binder plot, Uj{L) (Eq. 11), is given for several system sizes L. We estimate the critical 
temperatnre of the ferromagnetic-like phase transition from the crossing point: = 0.238. 

The value of Uj{L) at the point agrees well with that of the mean-hied model: Uf {L) = 
1 — r"^(l/4)/247r^ = 0.27...^^^*^ Thus we conclude that the elastic model with J > Jq, which 
shows the ferromagnetic-like transition belongs to the mean-hied universality class. This is 


consistent with a previous study.^^ (In Sec. IV, we discuss in detail the situation when J is 
closer to Jq.) 


B. Antiferromagnetic-like phase transition 

Next we study the case of the antiferromagnetic-Iike phase transition. For J < —0.02 
we expect that an antiferromagnetic-like transition takes place. Here we study the case of 
J = —0.1. Snapshots of the conhguration for L = 50 are given in Fig.|^(a). In contrast to the 
case of ferromagnetic-like order, the clusters consist of alternating LS and HS conhgurations. 
To hnd the clusters more easily, we also plot masked conhgurations of Fig. (a) in Fig. 
(b). Here the masked conhguration is given by cr' = (—1)®^+®!' x cTj and black (green) circles 
denote cr' = 1 (—1). There, small antiferromagnetic-like domains are observed at T = 0.29 
(left) and large antiferromagnetic-like domains at the middle (T = 0.25 which is close to 
the critical temperature T^. = 0.243). An ordered antiferromagnetic-like phase is observed 
at T = 0.2 (right). 

The Binder plot is given in Fig. (c) for several system sizes L. We hnd Tc = 0.243 from 
the crossing point, and the value of U^^{L) at the point agrees well with that of the short- 
range Ising model: Uj{L) = 0.61...^ Unlike the case of the ferromagnetic-like transition, 
the elastic model with J < Jq showing the antiferromagnetic-like transition belongs to the 
short-range Ising universality class. 


C. Correlation function 

We investigate the correlation function C{i,j)/C{l, 1), in which the value C(l, 1) is chosen 
as the reference value, and study the temperature dependence of the ordering patterns near 
Tc in both cases of the ferromagnetic-like and antiferromagnetic-like transitions. 

Figures]^ and [^illustrate prohles of the correlation function when J = 0.01 and J = —0.1, 


10 






respectively. The former shows the change of the correlation function in ferromagnetic-like 
transition, while the latter shows it in antiferromagnetic-like transition. In both transitions, 
(a) depicts a configuration below Tc, (b) a configuration close to Tc, (c) and (d) configurations 
above Tc. 

The conhgurations (c) and (d) in Fig. [^demonstrate that the correlation is still large at 
long distances in the ferromagnetic-like case, which is a characteristic of long-range interac¬ 
tion systems.^ In contrast to this feature, the correlation decreases rapidly at long distances 
in the antiferromagnetic-like case as shown in Figs [^ (c) and (d). This is characteristic of 
short-range interaction models like the Ising model. This analysis of the correlation func¬ 
tions is consistent with the analysis of the Binder cumulants, i.e., the universality class is 
different in the two cases. The difference between the mean-field and Ising universality 
classes is conhrmed by the correlation function. 

What causes the difference of the universality class between ferromagnetic-like and 
antiferromagnetic-like ordering? If we consider a conhguration in which the LS phase and 
HS phase coexist with a domain wall, the interface between the two phases causes an energy 
cost of the order of O(L^), which is the same mechanism as L^-dependence of the activation 
energy in the macroscopic nucleation of elastic interaction systems in 2D.^ (In the case of 
3D, the energy cost is of order 0{L^).) In this case, distortions of the lattice are propor¬ 
tional to the size of the cluster. Therefore, ferromagnetic-like large clustering is impossible 
because of the huge energy cost. As a result, the uniform configuration during the transition 
is universal in ferromagnetic-like ordering. 

The situation is different in the antiferromagnetic-like case. In appearance of the 
antiferromagnetic-like phase, the symmetry is broken between one configuration of LS, HS, 

LS, HS.. and the other conhguration HS, LS, HS, LS. If we consider a joint system of 

these two antiferromagnetic-like phases, the interface energy is of the order of 0{L) because 
these two antiferromagnetic-like phases have the same unit area in 2D (In the case 3D, the 
energy cost is of O(L^) because of the same unit volume). Thus, the size difference of the LS 
and HS molecules essentially does not cause energy cost. Therefore, antiferromagnetic-like 
ordering accompanied by clustering is possible as is the case of usual phase transitions in 
short-range interaction models. 
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FIG. 3: (color online) (a) Snapshots of the configuration for J = —0.1 and L = 50. Small clusters of 
the antiferromagnetic-like phase appear at the left (T = 0.29), large clusters appear at the middle 
(T = 0.25) near the critical point Tc = 0.243, and antiferromagnetic-like ordered phase is realized at 
the right (T = 0.2). (b) Snapshots of the masked configuration of (a), (c) Temperature dependence 
of for various system sizes L when J = —0.1. The horizontal dotted line = 0.61) denotes 
the fixed-point value of the Binder cumulant for the Ising universality class. 
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(C) 


(d) 


FIG. 4: (color online) Profiles of the correlation function of the system with J = 0.01 and L = 20 
around Tc = 0.238. (a) T = 0.200, (b) T = 0.240, (c) T = 0.260 and (d) T = 0.270. The diameter 
of each disk corresponds to the value of C{i, j)/C{1, 1). Yellow (Magenta) means the plus (minus) 
sign of the correlation function. 
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FIG. 5: (color online) Profiles of the correlation function C{i,j)/C{l,l) of the system with J = 
-0.1 and L = 20 around = 0.243. (a) T = 0.205, (b) T = 0.245, (c) T = 0.265 and (d) T = 0.275. 
The diameter of each disk corresponds to the value of C{i, j)/C{1,1). Yellow (Magenta) means 
the plus (minus) sign of the correlation function. 
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FIG. 6: Phase diagram in terms of the short-range interaction J vs temperature T. F denotes the 
region of the ferromagnetic-like phase, AF the region of the antiferromagnetic-like phase, and D 
the region of the disordered phase. The point J = Jq = —0.02 at T = 0 (downward-triangle) is 
the critical point between the ferromagnetic-like and antiferromagnetic-like phases at T = 0. Solid 
circles and squares denote the critical points for the antiferromagnetic-like and ferromagnetic-like 
transitions, respectively. The dashed line shows the critical temperature given by only the Ising 
interaction J. Open squares and circles denote the end points of the metastable ferromagnetic-like 
and antiferromagnetic-like phases (see Fig. [^, respectively. 

IV. PHASE DIAGRAM 

We depict the phase diagram J vs. T in Fig. Solid circles and solid squares denote 
critical temperatures (Tc) of antiferromagnetic-like and ferromagnetic-like phases, respec¬ 
tively. As pointed out in sec. [IT], the original point between the antiferromagnetic-like and 
ferromagnetic-like phases in the ground state is J = —0.02. It is worth nothing that the 
antiferromagnetic-like phase transition occurs at the point J = Jq{= —0.02) although the 
ferromagnetic-like and antiferromagnetic-like phases are degenerate at T = 0 for Jq. We 
hnd that the antiferromagnetic-like phase transition occurs even for J < —0.017, larger than 
J = Jq. It is considered as one of the reasons that the average density of the disordered 
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phase is closer to that of the antiferromagnetic-like phase than that of the ferromagnetic-like 
phase, and thus the generation of antiferromagnetic-like clusters is easier to realize than the 
generation of the ferromagnetic-like phase which is accompanied by global volume change. 

Considering that the critical temperature of the pure Ising modeP^ ( Eq.Q) (for both 
ferro and antiferromagnetic cases) is given by 


Tr = 


ln(l + V2) 


|J|^ 2.271 J|, 


(13) 


it is found that the critical temperatures of the antiferromagnetic-like transition, shown by 
solid circles in Fig. are a little bit larger than the temperatures T = 2.27| J| which is 
shown by the dashed lines in Fig. The critical temperature of the antiferromagnetic-like 
transition can be approximated by 


2.27(1 Jl+C), 


(14) 


where C ~ 0.01. Here it should be noted that Jo is just the dividing point between the 
ferromagnetic-like and antiferromagnetic-like phases at T = 0, and the critical temperature 
Tc is a function of J (not J — Jq) because the interface energy is a function of J. 

We consider the origin of the shift C of the critical temperature. We may attribute it 
to the cost of the elastic interaction at the interface. At the interface in a; or y direction 
between two phases (LS, HS, LS, HS... and HS, LS, HS, LS, ...), LS-LS molecular pair of 
the nearest neighbors (call LS-LS bond) and HS-HS molecular pair of the nearest neighbors 
(HS-HS bond) align alternately. The ideal intermolecular distance which gives the min¬ 
imum energy is different for LS-LS and HS-HS pairs. The surface energy due to the elastic 
interaction is calculated in a simple approximation as follows. HS-HS pair favors the dis¬ 
tance Tij = 2R}i but LS-LS pair favors Vij = 2i?L- Assuming that rtj = i?L + Rn is realized 
as a result of compromise, the energy costs per pair is AE ~ ^[rij — (i?j(cTj) -|- Rj{o'j))]‘^ ~ 
^[(i?L + i?H) -2i?L]" = ^[{RL + Rii)-2Rnf = !f[Rn-RL?=0.2, which leads to C ~ 0.1. 
However, the true value of C is much smaller and the relaxation of the conhguration would 
be necessary. 

Then we estimate the elastic interface energy as follows. We calculated the differ¬ 
ence of the elastic energies (Eq. between the system with no interface (complete 
antiferromagnetic-like phase) and that with an interface. For L = 20 and L = 40, we 
obtained relaxed conhgurations for both systems at T = 0 and estimated the elastic en¬ 
ergies. Dividing the difference of the elastic energies by the number of LS-LS and HS-HS 
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bonds on the interface, we have 6E ~ 0.023 (per interface bond) in the case of the interface 
in X (or y) direction for both L = 20 and L = 40. We also estimated 6E in the case of the 
interface in the diagonal direction, where interface line consists of LS-LS bonds or HS-HS 
bonds, and found 6E much smaller. Thus C is considered the order of 0.01. The energy 
cost due to interfaces is released to elastic interactions around the interfaces. The elastic 
interactions enhance antiferromagnetic-like ordering weakly. 

In contrast, in the ferromagnetic-like transition, the effective long-range interaction sig- 
nihcantly enhances the transition temperature and the values of are much larger than 
those of the antiferromagnetic-like transition. Unlike the antiferromagnetic-like transition, 
the critical temperature increases much more steeply than 2.27J with the increase of the 
value of J. This indicates that the elastic interaction and the short-range interaction J 
enhance ferromagnetic-like ordering synergetically with non-linear dependence. 

It is expected that around Jq the ferromagnetic-like and antiferromagnetic-like orders are 
nearly degenerate which causes a metastable structure of the ordered states. We study the 
metastable regions by the analysis of the dynamics of the order parameters under a sweep 
of J (Figs. (a) and (b)). 

In Fig. (a), we set ferromagnetic-like phase as the initial state at temperatures (T = 
0.01,0.03), and decreased gradually the value of J from —0.01 to —0.07, and observed the 
relaxation of {m?). We identihed the point where {m?) decreased large as the end-point of 
the ferromagnetic-like metastable region. Here, in the process of decreasing of the value of 
J, we changed J in steps of 0.001 and used 1000,000 MCS to equilibrate and the following 
1000,000 MCS to measure {m?) for each J. We determined the boundary of the metastable 
region where {rn?) becomes smaller than 0.9. We show the data of average over hve trials 
with the use of different random number sequences for L = 20 and L = 40. We do not hnd 
strong dependence on L for the metastability and conclude that the border of metastable 
state is well dehned. These points are plotted by open squares in Fig. 

In Fig. (b), we set antiferromagnetic-like phase as the initial state at temperatures 
(T = 0.01,0.03), and increased gradually the value of J from —0.03 to 0.03, and observed 
the relaxation of We identihed the metastable antiferromagnetic-like phase in the 

same way for the ferromagnetic-like case. The estimated boundary points are plotted by 
open circles in Fig. We hnd that metastability is realized and the region expands at lower 
temperatures. 
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FIG. 7: (color online) (a) Metastability of the ferromagnetic-like phase for the values of J, (b) 
Metastability of the antiferromagnetic-like phase for the values of J. + and x denote the order 
parameters at T = 0.03 and T = 0.01, respectively for L = 20. V Q denote the order 
parameters at T = 0.03 and T = 0.01, respectively for L = 40. 

V. OUTLINE OF SPIN CROSSOVER TRANSITIONS 


To study the characteristics of two-step (HS t-)- AF ^ LS) SC transitions on the basis 
of the analyses in the previous sections, we briefly summarize in this section the essential 
features of the types of SC transitions making use of the phase diagrams for the conventional 
(Ising-like) models. We will discuss the difference of the critical properties between the elastic 
model and Ising-like models in Sec. |VI[ 

Ising-like models have been developed to describe various SC phenomena. The Ising 
model with multi-fold degeneracy was studied for single-step (LS -H- HS) spin-crossover 
transitions to catch the entropy induced transitions.^ The model is given by 



Here Jp > 0 denotes the interaction (not magnetic) between the nearest neighbor molecules 
and D > 0 denotes the energy difference between the LS state (cTj = — 1) and the HS state 
(cTj = 1). The numbers of degenerate states u and r are associated with the LS and HS 


molecular states. The Hamiltonian (15) is equivalent to the following one in the partition 
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function.^ 


1 

'^ + 2 


iai + -{D - ksTIng)'^ai, ai = -1, 1, 


(16) 


T~i — —t/p^^CTjC 

{id) 

where T is temperature and g = ^ > Q. The second term is l-Les in Eq. ([^. This is the Ising 
model with an effective ffeld h{T) = — | (H — k^^Tlng). 

Using a phase diagram (T, h) for the Ising model shown in the upper panel of Fig. 
(a), we can obtain a better outlook to discuss the feature of the single-step SC transition. 


Raising the temperature in the model (15) causes a change of the order parameter, (m) = 
^ along a line h{T) in the phase diagram of the Ising model. The high spin fraction 

r (m)+l 

IS /hs = 

We define as the critical temperature of the Ising model without the field. In the mean- 
field theory, = zpJp, where zp is the coordination number. We also define Tcross(= kjing '^ 
as the temperature at which the effective field vanishes, i.e., h{T) = 0. Depending on the 
relation between and Tcmss, the transition between the LS and HS phases is classified into 
three cases.^ Typical three cases of temperature dependence of (m) are shown in the lower 
panel and corresponding paths (oblique lines) of h{T) are given in the upper panel. We 
find (I) gradual change when T^. < T^ross and (II) discontinuous (first-order) transition when 
Tc > Tcross- The second-order transition (case III) is realized only when = T^mss, which 
corresponds to the middle case in the lower panel. Considering the relation between h{T) 
and the spinodal lines (blue lines in the upper panel), we can classify the types of transitions 
in more detail.^ 

Ising-like models have been also extended to study two-step SC transitions and several 
important aspects were successfully clarified.®®^ The Ising-like model for two-step SC in 
bipartite lattices is given by 


H = - Tf ^ afaf -Jf ^ afaf 

(ieAjsA) (ieBjeB) 

- jaf ~ 

(is A, j SB) i 

where h{T) = —^{D — k^^T In g). Here A and B denote equivalent sublattices and Jf 
is the ferromagnetic-like intra-sublattice interaction and Jaf is the antiferromagnetic-like 
(Jaf < 0) inter-sublattice interaction. 

First we review the case of pure antiferromagnetic interaction ( Jf = 0). We depict a phase 
diagram (T, h) of the antiferromagnetic Ising model by the mean-field theory in the upper 
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FIG. 8 : (color online) (a) Phase diagram of the Ising model for single-step SC transitions (upper 
panel) by the mean-field theory and three typical cases of temperature dependence of (m) (lower 
panel). Spinodal curves are shown by blue lines and coexistence line by a green line. The critical 
temperature Tc is given by Tc = zpJp in the mean-held theory. Three cases of hrst-order, second- 
order transition and gradual change of (m) are shown in the lower panel. The corresponding path 
of h{T) for each case is also shown in the phase diagram (upper panel), (b) Phase diagram of the 
antiferromagnetic Ising model for two-step SC transitions and typical temperature dependences of 
(m) (lower panel). The red line denotes the critical line for the antiferromagnetic order {m\F), in 
which second-order transitions occur except T = 0. The dotted line gives (ttia) = 0 and (ttib) 7 ^ 0 
((me) = 0 and (ttt-a) 7 ^ 0). The critical temperature Tc at h = 0 is given by Tc = 2 :af|<A.f| in tbe 
mean-held theory. Temperature dependences of (m) for a two-step continuous transition (HS to 
AF to LS) and a one-step continuous transition (HS to AF) are drawn in the lower panel and the 
corresponding paths of h{T) are given in the upper panel. 

panel of Fig. [^(b).I^Two intrinsic order parameters (rriA) and (me), associated with the two 
sublattices A and B, exist, i.e., (mA) = (^b) = correspond 

respectively to the net magnetization per site in the sublattices A and B, respectively. N\ 
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FIG. 9: (color online) (a) Phase diagram of the Ising model with antiferromagnetic and ferromag¬ 
netic interactions for two-step SC transitions (upper panel) and typical temperature dependence of 
(m) (lower panel) of a two-step SC transition with first-order (lower T) and second-order (higher 
T) transitions. The critical temperature Tc at /i = 0 is given by Tc = zaf|>.^af| + zpJp in the 
mean-field theory. Tricritical points are shown by Pt{Tt,ht) and Qt{Tt, —ht). The blue lines de¬ 
note the limit of the metastability of the antiferromagnetic-like phase and the upper (lower) red 
line for 0 < T < Ft denotes the limit of the metastability of the HS (LS) phase. The red line for 
Tt < T < Tc is the critical line for antiferromagnetic-like order (ruAF)- (b) A two-step SC tran¬ 
sition with double first-order transitions, (c) Another type of two-step SC transition with double 
first-order transitions. Unlike the case (b), when Jp is relatively large, the decay temperature of 
the metastable HS phase can be lower than that of the metastable LS phase. 

{Nb) is the number of the sites in A (B) sublattice and = N^- We dehne two order 
parameters as (m) = ^"^a)+(™b) (mAn) = . The former and latter correspond 

to ferromagnetic-like and antiferromagnetic-like orders. 

The red line denotes the critical line in which the antiferromagnetic-like (AF) order (mAF) 
appears and it is the border between the regions of (ttia) = (ms) and {rriA) 7 ^ which 
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causes second-order phase transitions. Only at T = 0, the transition is discontinuous. The 
dotted line gives {mx) = 0 and {m-o) ^ 0 ((uib) = 0 and (ruA) ^ 0). Thus, only continuous 
(second-order) SC transitions are realized. The lower panel of Fig. (b) shows examples of 
a two-step (HS to AF to LS) continuous transition and a one-step (HS to AF) continuous 
transition, and corresponding paths (oblique lines) of h{T) are also drawn in the upper 
panel. 

If Jp 7 ^ 0, the metastable region of the antiferromagnetic-like phase and that of the 
ferromagnetic-like phase appear at low temperatures and they exist in the region 0 < T < Tt, 
where Tt is the temperature of tricritical points Pt{Tt,ht) and QtiTt, —ht) (see the upper 
panel of (a)). For 0 < T < Tt, the upper and lower blue lines denote the limit of 
the metastability of AF phase (i.e., AF spinodal line), while the upper (lower) red line 
corresponds to the limit of the metastable HS (LS) phase (i.e., HS (LS) spinodal line). The 
limit of the held for the metastable AF phase at T = 0 is h = ^af|A 4 f| + zpJp = he 
and h = —he, and the limit of the held for the metastable HS (LS) phase at T = 0 is 
h = -Saf|A 4 f| — ^f-^f = ho {h = —ho). The critical temperature (Neel point) Tc at h = 0 is 
given by Te = zafI-^afI + ^^f-^f in the mean-held theory. Here zaf (^^f) is the coordination 
number for the inter(intra)-sublattice nearest-neighbor sites. In the region Tt < T < Te, the 
red line denotes the critical line of the antiferromagnetic-like order tuaf for second-order 
phase transitions as the same as the case Jp = 0. 

The lower panel of Fig. (a) illustrates an example of a two-step transition with hrst- 
order (lower T) and second-order (higher T) phase transitions. In Fig|^(b), an example of a 
two-step transition with double hrst-order transitions is depicted. Depending on the relation 
between two tricritical points {Pt and Qt) and the line h{T), the two-step transitions are 
classihed into three cases: (I) hrst-order and second-order, (H) double hrst-order, and (HI) 
double second order. When the line h(T) locates above the two tricritical points, the type 
(H) is realized, while h{T) locates below the two tricritical points, the type (HI) is realized. 
When h{T) locates between the two tricritical points, the type (I) occurs. Figsj^ (a) and 
(b) correspond to the types (I) and (H), respectively. 

When Jp is so large as 2 ;pJp > 2 ;af|</af|, another pattern of two-step transitions of 
double hrst-order transitions can be realized as depicted in Figj^ (c). Here the limit of the 
metastability of the HS (LS) phase (red line) appears at the lower (h < 0) (upper (h > 0)) 
half plane, and thus the decay temperature of the metastable HS phase can be lower than 
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the decay temperature of the metastable LS phase. In SC compounds, this case of double 
hrst-order transitions has not yet reported in experiments but it may be found in the future. 
The three patterns of two-step transitions were also studied by Monte Carlo methods.SH In 
the elastic interaction model (Eq. [^, the elastic interactions {ki and ^ 2 ) and the short-range 
interaction J take the roles of the interactions Jp and Jaf of the Ising-like models. 

VI. SUMMARY AND DISCUSSION 

We obtained a phase diagram as a function of the short-range interaction J, where we 
found that both the ferromagnetic-like and antiferromagnetic-like transition temperatures 
are enlarged by the elastic interaction. However, we found that the natures of ferromagnetic¬ 
like and antiferromagnetic-like phase transitions are qualitatively different. In the case of 
the ferromagnetic-like phase transition, the elastic interaction acts as an effective long-range 
interaction and it signihcantly enhances the ferromagnetic-like ordering, where the system 
belongs to the mean-held universality class. The increase of the critical temperature is 
much larger than that expected from the Ising interaction. The synergetic effect of the 
elastic interaction and the short-range interaction amplihes the ferromagnetic-like ordering 
and causes high critical temperatures. 

In sharp contrast to this case, in the case of the antiferromagnetic-like phase transition, 
the system belongs to the Ising universality class. In this case the long-rang interaction 
due to the elastic interaction is irrelevant, and clustering of the ordered phases is observed. 
The elastic interaction raises the critical temperature of antiferromagnetic-like order, as 
well. This is because the staggered structure is stable for the elastic interaction, but the 
contribution to enhancement is small. The antiferromagnetic order is mainly enhanced by 
the short-range interaction. 

We also conhrmed different natures in the conhgurations of the correlation function near 
the critical point in the cases of ferromagnetic-like and antiferromagnetic-like phase tran¬ 
sitions. Besides, we found that metastable ferromagnetic-like and antiferromagnetic-like 
phases exist near the region in which both orders are nearly degenerate. 

The present study was performed for the two-dimensional model, but the conclusion can 
be extended to cases of three dimensions because the physical mechanisms studied in this 
paper is the same as in three dimensions. 
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In the present study, we focused on the critical properties of the model along the co¬ 
existence line, i.e., Hefr = 0. In the SC system, Heff changes with the temperature, and 
reflecting the structure around the coexistence line, the system can exhibit a two-step SC 


transition with temperature change. We depict two types of two-step transitions in Figs. 10 


(a) and (b). The former shows double continuous transitions and the latter shows hrst-order 
(lower T) and continuous (higher T) transitions. In this Monte Carlo simulation, we used 
1000,000 MCS to equilibrate and the following 1000,000 MCS to measure (m). From the 
results of this study, it is found that the SC phase transition between the HS phase and the 
intermediate phase (antiferromagnetic-like phase) can be a second order phase transition of 
the short-range Ising interaction type, while the SC phase transition between the HS and LS 
phases is of the mean-held type when it is of the second order. The approach of the present 
study enables us to capture new aspects of two-step SC transitions from the view point of 
the elastic interaction and the short-range interaction. 

The phase diagrams of HS, LS, and AF phases for the elastic interaction model with the 
short-range interaction may have similar features to those obtained by Ising-like models in 
the mean-held theory in Sec. |V] (see Fig. and Fig. |^. However, the nature of the phase 
transition of the model is not the same as that of the Ising-like models. The universality 
classes of the transitions are diherent from those of the Ising-like models. The regions of 
the metastability for the HS, LS, and AF phases are determined by local interactions in 
the Ising-like models, while those are determined by both local and global stability in the 
elastic model. Thus the quantitative features of the phase diagrams in the elastic model are 
different from those in the Ising-like models, which is non-trivial. The details of the phase 
diagrams and correspondence to SC compounds will be studied in the future.!^ 

When J = Jq, the antiferromagnetic-like phase appears although the ground states 
(T = 0) of the two phases are degenerate at this value of J. This fact indicates that the 
antiferromagnetic-like phase is favorable at hnite temperatures. In the text we considered 
a reason why the antiferromagnetic-like phase is realized easier than the ferromagnetic- 
like phase. We may attribute it to a dynamical effect. The average fraction of HS and 
LS molecules of this ordered phase is almost the same as in the disordered phase and the 
antiferromagnetic-like phase can be realized without modihcation of the density. We may 
also consider this fact from the view point of entropy, i.e., a kind of mechanism of order by 
disorder. The detailed study will be done in the future. 
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(b) 


FIG. 10: (color online) (a) Two-step SC transition with double continuous (LS -H- AF and AF <-)• 
HS) transitions. D = 1.6 and g = 1000. The symbols Q A denote (m) for L = 20 and L = 40, 
respectively, (b) Two-step SC transition with first-order (LS -H- AF) and continuous (AF -H- HS) 
transitions. D = 1.02 and g = 1000. 
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FIG. 11: (color online) The definition of Ojn. 


Appendix A: Interactions to maintain the square symmetry and the coexistence of 
high-spin and low-spin phases 

If only the nearest neighbor interaction (Eq. ([^) is used, the lattice may distort into 
a rhombic shape. To prohibit such deformation, we need some additional interaction. In 
the text, we applied the next nearest interaction "Hnnii (Eq. (|^). We may consider other 
interactions instead of the choice of A kind of interactions to maintain the angle 

between bonds for desired lattices has been frequently used.^^^® Thus in the same way we 


focus on the angle (Oju) between the bonds i — j and i — I (Fig. 11), which is dehned by the 
relation 


cos 9jii = 




(Al) 


To maintain a square lattice, i.e., Oju equal to 7r/2, the following interaction can be adopted. 




(A2) 


where C* is a positive constant, and the summation is taken over all pairs of bonds. 

For the conhguration of the ferromagnetic-like phase or antiferromagnetic-like phase in 
Fig-B l-Le takes the minimum value 'He = 0, and this term plays a role in maintaining 
the square symmetry. In this interaction, the ferromagnetic-like and antiferromagnetic-like 
phases have the same energies at T = 0 and both phases are the ground state. The origin 
of J is zero (Jq = 0) in this case. 

Unlike the case of Eq. (|^, in this choice (Eq. (A2)), the conhguration entropy is different 
for HS molecules and LS molecules, as analyzed below. We consider the motion of molecules 
around the position in the complete LS (HS) phase (as in Figs (a) and (b)). Let the 
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fluctuation of a molecule in a:-direction (?/-direction) defined as dx {dy). The potential term 
for the molecule U, given by cos^ dju, has the form 


U oc 


{dx + dyY 

w 


(A3) 


in the leading term of dx and dy, where R is the molecular radius. Because R is larger in the 
HS state and the entropy gain is larger in the HS state, the HS state is more favorable than 
the LS state. Indeed, we observed that only HS state is realized in the simulation. LS and HS 
states are not symmetric anymore. This situation is similar to the Ising model with nonzero 
magnetic field. Thus critical phenomena do not occur in this case. A field for cancelation 
of this difference was artificially applied when critical properties were studied.^^^^ Thus, 
we adopt the present choice to avoid this complication, where Rnnn has no /^-dependence 
as dx'^ -|- dy'^ in the leading term and the system has a critical point. Rigorously speaking, 
the HS state is more favorable than the LS state even in this case. However, Rnnn has 
R-dependence at higher order and the difference is small and can be ignored in the practical 
calculation (no influence on simulations). This analysis holds in three dimensions. Namely, 
Eq. (A3) has the dependence of ^ or for the fluctuation {dx, dy, dz), 

while Eq. ([^ has no i?-dependence as dx'^ -f- dy'^, dy‘^ + dz'^, or dz'^ + dx"^. 

In order to confirm the above mentioned effects, we performed Monte Carlo simulations 


with the interaction (A2) with (7 = 4. For J < 0 we found that antiferromagnetic-like 
transition occurs and it has the same critical properties as in the case of the text. Namely, 
curves of T dependence of cross for different system sizes L and the value of at 
the crossing point is = 0.61 (not shown). For example, is estimated = 0.138 


when J = —0.04. It is clear that the interaction (A2) has the same contribution to the two 
coexistent anitiferromagnetic-like phases and has no effect on the critical properties. 

In contrast, for the ferromagnetic-like region, it has a big effect on the critical property. In 


Fig. 12 (a) the temperature dependence of {m?) is given for J = 0.04. Here we find a sharp 
change of magnetization {m?) which indicates ferromagnetic-like phase transitions, and the 
magnetization seems to have a critical point. However, if we plot the Binder cumulant for 
the ferromagnetic-like transition, a strange dependence on the size was found. In contrast 
to the antiferromagnetic-like case, curves of T vs. for different system sizes do not cross 


in this ferromagnetic-like case, shown in Fig. 12 (b). Thus we conclude that the critical 
point does not exist in this case. When J = 0.001, an antiferromagnetic-like transition 
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(a) 


(b) 


FIG. 12: (color online) Temperature dependence of {m?) (a) and temperature dependence of Uj 
(b) for various system sizes L when J = 0.04. 


takes place, and for J > 0.01 the system shows a ferromagnetic-like transition. There is 
a critical valne of J (Jc) between the ferromagnetic-like and antiferromagnetic-like phases. 
This critical value locates in the region of 0.001 < Jc < 0.01. This is considered a similar 
hnite temperature effect which is discussed in the text. 

As another choice, we may adopt a next-nearest neighbor interaction to realize Jq = 0. 

Hnnn = ^ Y. ■ (A4) 

((i,k)) 

Here the next nearest neighbor interaction i — k depends on the bonds i — j, j — k, k — I 
and I — i, and we dehne 


‘^Rijkl — '[{Ri{(^i) + Rji^^j)) + {Rji^j) + Rk{(^k)) + {Rk{(^k) + Rl{(^l)) (AS) 
+ {Ri{ai) + R,{ai))]/4 
= [i?j((Tj) -|- Rjicj) -|- Rk{ak) + Ri{(Ji)\/2. 


In this case the states of Fig. [^(ferromagnetic-like and antiferromagnetic-like conhgurations) 
are the ground state and Jo = 0. However, this model gives a similar situation of non-crossing 
of Uj-T curves although some cases of L show crossing. This reason is not the same as the 


case (A2). The reason is not so clear but 4-body interactions may cause such complex 


dependences. Thus, although Jq ^ 0, we adopted Eq. (j^ for the purpose of maintaining 
the square lattice in the present work. 
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